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Abstract. The orthogonal projection of a 4-cube onto a uniform random 
3-subspace in M 4 is a convex 3-polyhedron P with 14 vertices almost surely. 
Three numerical characteristics of P - volume, surface area and mean width - 
are studied. These quantities, along with the Euler characteristic, form a basis 
of the space of all additive continuous measures that are invariant under rigid 
motions in M 3 . While computing statistics of {vl, ar, mw}, we encounter the 
generalized hyper geometric function, elliptic integrals and Catalan's constant. 
A new constant 7.1185587167... also arises and deserves further attention. 

A planar shadow of a 3-cube {4-cube} is a convex hexagon {octagon} almost 
surely. In an earlier paper [lj, joint moments of hexagonal {octagonal} area and 
perimeter were computed. It is natural to speculate on n-cubes, under the action of 
corank-1 projections rather than rank-2. 

Let C denote a 4-cube with edges of unit length, centered at the origin. To 
generate a random 3-subspace S in M 4 , we select a random point U uniformly on the 
3-sphere of unit radius. The desired subspace is the set of all vectors orthogonal to 
U. 

We then project the (fixed) 4-cube C orthogonally onto S. This is done by forming 
the convex hull of images of all vertices of C. The resultant polyhedron in the 
hyperplane has 14 vertices almost surely. More precisely, if U — (x, y, z, w) is of unit 
length, then the matrix 



Ma 



/ , /l _ ~2 X V xz xw \ 
z 2 +w 2 y_z yw_ 

1 ~ X2 y/(l-X 2 )(z 2 +W 2 ) y/(l-X 2 )(z 2 +W 2 ) 







ir 



y/z 2 +w 2 \Jz 2 +w 2 

\ J 



projects C orthogonally onto a hyperplane, rotated in R 4 to coincide with the 3- 
subspace spanned by (1, 0, 0, 0), (0, 1, 0, 0) and (0, 0, 1, 0) for convenience. Let T be 
the 4-row matrix whose columns constitute all vertices of C. Then the first 3 rows 
of M4T constitute all images of the vertices in M 3 and 3-dimensional convex hull 
algorithms apply naturally. Such calculations provide the underpinning for our work. 
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Spherical coordinates in M 4 : 

x = cos 8 simp sin ip , y = sin8 simpsinip, z = cos<£> simp, w = cosip 

will be used throughout for U, where < 9 < 2n, 0<ip<7T,0<ilj<7i. The 
corresponding Jacobian determinant is sin ip sin 2 ip; it is best to think of (8, ip, ip) as 
possessing joint density ^ sin ip sin 2 ip. 

1. Volume 
The projected polyhedral volume vl in R 3 is equal to 

\x\ + \y\ + \z\ + \w\ , 

given a unit vector U = (x,y,z,w). It follows that 

tt/2 it/2 tt/2 

E (vl) = 4 ■ 16 J J J cosip ^— j sin cp sin 2 ip dip dtp dO 
ooo 

= — = 1.697652726313550... 
3tt 



tt/2 tt/2 7r/2 

E (f/ 2 ) = 4 • 16 J J J cos 2 ip — - sin y9 sin 2 ip dip dip d9 
ooo 

tt/2 tt/2 tt/2 

+ 12-16 J J J cos ip sin ip cosip sin ip sin 2 ip dip dtp d9 
ooo 

= 1 + - = 2.909859317102744.... 

7T 

More generally, starting with a unit n-cube, the projected (n — l)-polyhedral volume 
vl in M n_1 satisfies 

n r(f) 



E(vl) 



vrr(H±i)' 
2(n - 1' 



E (vl 2 ) = 1 + 

These volume moment formulas are the same as formulas in [2] for the mean width 
and mean square width, respectively, of the n-cube itself. Such duality is a special 
case of a theorem proved in [2HHE1E]- We also have max (vl) = y/n and min (vl) = 1. 
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2. Surface Area 
The projected polyhedral surface area ar in M 3 is equal to 

2 ^a/x 2 + y 2 + \/ x 2 + z 2 + V 'x 2 + u> 2 + \J y 2 + z 2 + \/y 2 + w 2 + V z 2 + w 2 ^ , 

given a unit vector [/ = (x,y,z,w). It follows that 

tt/2 tt/2 tt/2 

E (ar) = 6 • 16 ■ 2^ ^ ^ ^/ cos 2 y? sin 2 ^ + cos 2 ip ^—^ sin </? sin 2 ip dip dip d9 = 8 
ooo 

and, since x 2 + y 2 = sin 2 sin 2 ■?/>, 

tt/2 tt/2 tt/2 

E (ar 2 ) = 6 • 16 • A J J J (cos 2 ip sin 2 ip + cos 2 ^) — - sin </? sin 2 ip dip dip d9 



ooo 

tt/2 tt/2 tt/2 



+24 • 16 • A J J J sin ip sin ip \j sin 2 9 sin 2 </? sin 2 ip + cos 2 ^ — ^ sin 99 sin 2 ^ dip dip d6 
000 

tt/2 tt/2 tt/2 

+6 • 16 • 4 J J J sin ip sin ip \J cos 2 <p> sin 2 ^ + cos 2 ip — ^ sin ip sin 2 ^ a 1 -?/; dip d9 
000 

= 12 + 6C 4 + 3tt = 64.136130261087789.... 

Details on derivation of the mean square result, including a formula for the constant 
£4, will appear shortly. We also have max (ar) = 6y / 2 and min (ar) = 6. 

More generally, starting with a unit n-cube, the projected (n — l)-polyhedral 
surface area ar in M n_1 satisfies 

E(ar) = ^ {n - 1)n r(f) 



2 rp±!)' 

E (ar 2 ) = 4(n - 1) + („ - 2)(n - 1)C B + { - - 3)(w ~ 2)(w " % . 
For the special case n = 3, we have 

^- = a/x 2 + y 2 + Vrr 2 + z 2 + \/ y 2 + z 2 = \/\ — x 2 + a/1 — y 2 + Vl — z 2 = ^ 

given that x 2 +t/ 2 + ^ 2 = 1, which implies a closed-form expression for ( 3 (using results 
from an upcoming subsection "2D Analog"). Numerical evidence strongly suggests 
that ( n = ( 3 for all n > 4, but a rigorous proof is not known. 
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3. Mean Width 
The projected polyhedral mean width mw in M 3 is equal to 



- (V 1 - x 2 + ^l-y 2 + Vl-z 2 + Vl - w 2 ) , 



given a unit vector U = (x,y,z,w). It follows that 

tt/2 tt/2 tt/2 



K(mw) = 4 • 16 • —J J J a/ 1 — cos 2 ip — - sin <p sin 2 ^ dip d6 



ooo 



tt/2 tt/2 tt/2 



E (m-u; 2 ) = 4 • 16 • i y ^ ^ ( 1 - cos 2 *0) 7^ sin </? sin 2 ip dip dip d0 



ooo 

tt/2 tt/2 tt/2 



+12-16- —J J J \J 1 — cos 2 y? sin 2 -^a/ 1 — cos 2 ip — ^ sin ip sin 2 ^ dip dip d6 



= 3 (- + - + - J = 2.883026903647544... < 1 + - = E (vl 2 ) . 
\4 8 ir/ 7r v ' 

We also have max (mw) = a/3 and min (mw) = 3/2. 

More generally, starting with a unit n-cube, the projected (n — l)-polyhedral mean 
width mw in M n_1 satisfies 

E (mw) = _ / 2 ^. = E (uZ) 

but a general formula for E (mw 2 ) is unknown. We examine the cases n = 3 and 
n = 5 in the following subsections. 
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3.1. 2D Analog. As outlined in [T], 

2 



- (Vl - x 2 + y/l-y 2 + Vl 

IT \ 

given that x 2 + y 2 + z 2 = 1, hence E (row) = 3/2 and 



7T/2 7T/2 

2 ' ' 



E (mw 2 ) = 3 • 8 • J J (1 - cos 2 y>) — sin ip dip d9 



2 X2 





7T/2 7T/2 



+6 • 8 ■ ^— J j J \J 1 — cos 2 6 1 sin 2 ipy/l — cos 2 ^ — sin <iyj <i# 



= ^ [4 + 3tt 3 F 2 (-§, |, |; 1, 2; l)] = 2.253091059149751... < 1+- 

where 

r(6i)r(6 2 ) A r>! + k)r{a 2 + fc)r> 3 + k) z k 



s^lai,^, a 3 ; 6i, b 2 ; z) - —————— ^ 



r(a 1 )r(a 2 )r(a 3 ) ^ r(6i + A;)r(6 2 + k) k\ 

is the generalized hypergeometric function. As a byproduct, the constant ( 3 defined 
earlier is equal to 37r 3 F 2 (— |, |, |; 1, 2; l). We also have max (mw) = 2y/Q/ir and 
min (mw) = 4/tt. 
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3.2. 4D Analog. We here have 



4 5 

mw = — > a/1 — x 2 
3vr V i 



3=1 



given that Y^j=\ x2 j = 1; hence E (mw) = 15/8 and 



E (mw 2 ) = 32 • (^j (51 + 20 J) 



where 

tt/2 tt/2 tt/2 tt/2 



J J J J (l — cos 2 </? 3 ) sin </? 1 sin 2 </? 2 sin 3 </? 3 d(p 3 d(p 2 d(p 1 d9, 





tt/2tt/2 tt/2 tt/2 



J = J J J J \J 1 — cos 2 ip 2 sin 2 (/9 3 a/1 — cos 2 y? 3 ^— ^ siny9 1 sin 2 <p 2 sin 3 ^3 d(p 3 d(p 2 dp 1 d9. 
0000 

The expression for E (mw 2 ) simplifies to 



I . . - ,iy l 



81tt 4 



144tt 2 - 10r (i) 4 + 45tt 3 (8 3 F 2 (-1, ±, |; 1, 2; l) - 3 F 2 (-§, |, f ; 1, 3; l)) 



= 3.516040901689803... < 1+-. 

7T 



We also have max (mw) = 8a/5/ (3%) and min (mw) = 16/(37r). The fact that E (mw) 
at one level becomes min (mw) at the next level is interesting. 



Rank-3 Projections of a 4-Cube 



7 



4. Details 
With regard to E (ar 2 ), we first prove that 

tt/2 tt/2 tt/2 

J J J sin (p sin ipyj cos 2 (p sin 2 ip + cos 2 ip -^—^ sin </? sin 2 ip dip d(p d6 — ■ 



After rearranging and integrating out the integral becomes 

tt/2 tt/2 

^ J cos ^ sin 3 ^ J a/ cos 2 y9 tan 2 ip + 1 sin 2 </? d</? d0 
o o 

tt/2 

1 /" . 3 (-1 +tan 2 ^)^(i tan^) + (l + tan 2 ^)i^(i tan^) 



= -J cos^sin # 

o 

where i is the imaginary unit and 

tt/2 tt/2 

k(0 = [ 1 de, e(0 = [ Ji-e*Ho?de 



are complete elliptic integrals of the first and second kind. Let t = t&nip, then 

dt = sec 2 ip dip = (l + t 2 ) dip 



hence 



3 ./, / 1 \ 3/2 / x 1/2 



cos ip sin 3 ip 3 / 1 \ 



tan 2 -0 



d0 = cos ip simp dip — 1 dip 



5/2 / \ 1/2 



\i + t 2 ) V i + *V 



-dt 



(i+t 2 r 



and the integral becomes 

oo 

1 f (-l + t 2 )E(it) + (l + t 2 )K(it) 



12tt J (1 + t 2 ) 

o 



3 



12tt 



y -^tdt-2 / K -^tdt+ / y -^tdt 

(1+t 2 ) 2 y (i+t 2 ) 3 y (i+t 2 ) 2 



.o o 

7T 7T 7T 7T 



96 128 192 128 
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The other integral 

tt/2 tt/2 tt/2 



sin p sin ip \/ sin 2 6 sin 2 p> sin 2 ip + cos 2 ip — - sin (p sin 2 ip dip dp d6 



where 



ooo 

is harder. We rearrange it as 

tt/2 tt/2 tt/2 

J J cos ip sin 3 ip J \J sin 2 6 sin 2 p tan 2 ip + 1 sin 2 p d9 dp dip 



tt/2 tt/2 

= — - J cos ip sin 3 ip J E (i sin p tan ip) sin 2 </? 




f(iP) = (-2 + 4tan 2 ^)E^ V / ^i^) , 

<?(V>) = 2 (1 - 2tan 2 ^ + sec^) E (i y/^^j K (i y/^^j , 

h(ip) = (-l + tan 2 ^-sec^)K J™£?J . 
Let t = a/ (sec ip — l)/2, then 1 + 2t 2 = sec ^ and 

4tdt = sec ipta.ni/j dip = (l + 2t 2 ) ((l + 2t 2 ) 2 - l)^ # 

hence 

cos ^ sin 3 ip 3,-,,, /" 1 V f-, 1 \ 1/2 ,, 
-rf^ = cos ipsmip dip = 1 2 dtp 



tan 2 V V 1 + 2i2 / V (1 + 2t 2 ) 2 

1 \V 1 \ 1/2 1 



" (4t Hl + 2t 2 J I 1 (1 + 2^ ((l + 2t 2 ) 2 -l) 1/2 ^ 
(1 + 2t 2 ) 5 



Rank-3 Projections of a 4-Cube 



9 



The integral becomes 



3tt 2 



-2 + 4((l + 2t 2 ) 2 -l 



E{itf 



-tdt 



3tt 2 



3tt 2 



oo 

/ 



oo 

/ 



(l + 2t 2 ) 5 

1 - 2 ((1 + 2t 2 ) 2 - l) +(l + 2t 2 )] E(it)K(it) 



(1 + 2t 2 y 



tdt 



-1+ ((l + 2t 2 ) 2 -l) - (1 + 2t 2 



K(it) 2 



(l + 2t 2 r 



-tdt 



which simplifies to 



oo 

4 f (8t 4 + 8t 2 -l)E{i tf - 2 (4t 4 + 3t 2 - 1) E (it) K (i t) + (2t 4 + t 2 - 1) K (i tf 



3tt 2 



/ 



(2t 2 + 1)' 



tdt. 



Call this expression ^4/256. The constant C4 = 7.118558716719735... coincides with 
( 3 to high numerical precision, but algebraic confirmation remains open. 

5. Correlations 
The joint moment for volume and surface area is 

tt/2 tt/2 tt/2 

E(vl-ar) = 12 • 16 ■ 2J J J cos ip \J cos 2 p sin 2 ip + cos 2 ip 



000 



+ cos ip \J sin 2 9 sin 2 ip sin 2 ip + cos 2 ~tp sin 2 if) 



6(1 + -)= 13.639437268410976... 

7T 



2tt 2 



sin (p sin 2 ip dip dip d6 



which, together with earlier results, implies that the correlation between vl and ar is 
0.945.... Likewise, we have 

tt/2 tt/2 tt/2 

E(vl-mw) = 4-16- -J J J cosip\/l — cos 2 ip — - sin ip sin 2 ip dip dp d6 

000 

tt/2 tt/2 tt/2 



+12-16.1 



J J J cos p sin ip \J 1 — cos 2 ip — - sin p sin 2 ip dip dp d9 



000 



= - + - = 2.886619772367581..., 

4 7T 
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tt/2 tt/2 tt/2 



E (ar ■ mw) 



12 • 16 • 2 ■ 








cos 2 (p sin 2 ip + cos 2 ip \/\ — cos 2 ip 

1 



+ \J sin 2 9 sin 2 <^ sin 2 -0 + cos 2 (p sin 2 ^ \fl — cos 2 -0 
3 (5 + 2G) 9tt 



2tt 2 



sin sin ip dip dip d9 



+ — = 13.592597187518807.. 
4 



where 



fc=0 



(2A: + 1) 1 



is Catalan's constant [7]. The correlation between vl and is consequently 0.870... 
and the correlation between ar and mw is 0.973.... 

6. Related Work 
Let us return to an open issue in pQ surrounding rank-2 projections of a 4-cube. On 
the one hand, octagonal perimeter is equal to 



2 (Vl -P 2 -x 2 + v 7 ! - q 2 - y 2 + VT 



r 2 — z 2 + Vl — 



s 2 — u» 2 



given orthogonal unit vectors £/ = (x, y, z, w) and V = (p, q, r, s). This formula makes 
accurate calculation of the second moment of perimeter 28.495... more feasible. The 
maximum value of perimeter is 4 a/2; the minimum value is 4. 

On the other hand, octagonal area does not appear to have so simple a description. 

Let 



Q 

r 



cos k sin A 



/ -y \ 

X 

V * / 



+ sin k sin A 



w 

X 

\-y J 



+ cos A 



f -w\ 

—z 

y 

\ x ) 



where < k < 27r, < A < n. Let 



«1 


= r y + s y 


— qz — sz — qw + rw 




= ry — sy 


— qz — sz + qw + rw 


«3 


= ry + sy 


— qz + sz — qw — rw 




h=pq 


— ps + xy — xw, 




b 2 = pq 


— pr + x y — x z, 
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&3 = pq + ps + xy + xw, 
c = 2(1 -p 2 -x 2 ). 

Then, in a neighborhood of (9 , ip , i[) , k q , A ) = (1, 1, 1, 1, 1), octagonal area is equal 
to 

- \a x b 2 + a 2 (c -bi + b 2 ) + a 3 &i] ; 
c 

in a neighborhood of (6*o, v^o; ^oi K o, Ao) = (|, |, §, §, §), octagonal area is equal to 

- [ai& 2 - a 2 (6i - 62) + a 3 (c + 61)] ; 
c 

in a neighborhood of (0 Q , tp , ip , k , A ) = (f , §, |, |, |), octagonal area is equal to 

i [ai& 2 - a 2 (c + foi — 6 2 ) + «3&i] ; 

in a neighborhood of (#0, (p Q , ^o> ^0) = (| ; §> §> f j §); octagonal area is equal to 

1 

- - Q i( c ~ h) ~ a 2 (h - b 2 ) + a 3 &i ; 
c 

in a neighborhood of (#0, V^O' ^o> K o, Ao) = (|, |, |, |, |), octagonal area is equal to 

- [aib 2 - a 2 (h - b 2 ) - a 3 (c - bx)] ; 
c 

in a neighborhood of (9 , (p , ip , k , A ) = (§, 1, §, §, |), octagonal area is equal to 

^ [ai(c + b 2 ) - a 2 (b 1 - b 2 ) + a 3 &i] ; 

Other branches also exist, but for reasons of space, we stop here. The maximum 
value of octagonal area is 1 + a/2 [51 |6]; the minimum value is 1. 

An addendum to [8] now exists, with voluminous detail on the derivation of a 
relevant density function [9]. 

The present work is the fifth (and final) in a pentalogy that began with [TO] , 
continued with [2j [TTJ and then again with [TJ. 
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8. Addendum 
Here is a proof that ( 5 = £ 4 . For n — 5, we have 



ar = 2 y x2 j + x l 

l<j<&<5 



given that X^=i x j = 1- The hard integral for n = 4 in a preceding subsection 
"Details" here becomes 

tt/2 tt/2 tt/2 tt/2 

l l l l ^ (/?2 ^ 3 \l S ^ n2 ^ S ^ n2 ^ l S ^ n2 ^2 S ^ n2 ^3 C ° S2 ^ 3 


3 

x — — sin ip 1 sin 2 </? 2 sin 3 </? 3 d(p 3 dip 2 d(p 1 d9 

tt/2 tt/2 tt/2 tt/2 

= J J J cos 9^3 sin 4 ^* \J sin 2 9 sin 2 <p x sin 2 y2 2 tan 2 (p 3 + 1 



3 

x — — sin 2 (p 1 sin 3 <y2 2 d# d(p 1 d(p 2 d(p 3 

OTT 2 

tt/2 tt/2 tt/2 

= -^-^ J J cos </7g sin 4 _E (i sin ip^ sin (^> 2 tan ( -p 3 ) sin 2 sin 3 y? 2 d(p-^dip2d(p 3 





tt/2 tt/2 



where 



— J J sm^ 2 cos^sm ^ ^2^ 



/(y, V>) = (-2 + 4 sin 2 y tan 2 ^ E ^ ^ , 
g(ip,ip) = 2^1 — 2 sin 2 ip tan 2 ijj + \J 1 + sin 2 </? tan 2 ^ 



x 



£ ^ ^/ y/ 1+sin 2 ip tan 2 V>- 1 ^ ^ ^ ^/ ^l+sin 2 y?tan 2 ^ 



■?/>) = ^—1 + sin 2 ip tan 2 ip — \J\ + sin 2 </? tan 2 ijj^j K (i\J 



\J 1+sin 2 <ptan 2 ip—1 
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Let 

then 

hence 



and 



/ \ l+sin 2 a? tan 2 ib~l 

u=\J x 2 > v — sm(p 

1 + 2m 2 = a/1 + sin 2 ip tan 2 ip, y/1 — v 2 = cos ip 

2 , (1 + 2m 2 ) 2 -1 4m 2 (1+m 2 ) 
tan 2 ip = ± / = K — 

9l 4m 2 (1+m 2 ) 4m 2 (1 +u 2 ) +v 2 
sec 2 ip = l + ^ - = i 



2(1 + 2m 2 ) (4m) — = (sin 2 ip) (2 tan ip sec 2 if;) 

^ /4m 2 (1 + m 2 ) 4m 2 (1 - - :rr- r 2 



(4m) [4m 2 (1 + m 2 ) + v 2 } VTTv 2 



It follows that 



du _ [4m 2 (1 + m 2 



v 2 } VT+rf 



2(1 + 2u 2 )v 



dip 



hence 
hence 

sin <p cos ip sin ip 
tan 2 ip 



d(u, v) 



d(ip,ip) 



[4m 2 (l + m 2 ) + v 2 ] VT+^VT^ 2 

2(1 + 2m 2 )v 



dip dip = sin <p> cos 3 ip sin 2 ip dip dip 



Au 2 (1 + m 2 ) + v< 



2\ -3/2 



1 - 



4m 2 (1 + u 2 )+v 



2\ -1 



d<p dip 



Au 2 (1 + m 2 



[4m 2 (1 + m 2 ) + v 2 f 2 4m 2 (1 + m 2 ) + v 2 

d(<p,ip) 



dip dip 



[AU 2 (1 + M 2 ) + V 2 f 2 



[4m 2 (1 + m 2 ) +v 2 ] 7/2 



Au 2 (1 + m 2 ) 



d(u, v) 



du dv 



a 2 P\ o 2(1 + 2m 2 )*; j 

4m 2 VI + m 2 v - ' -du dv. 
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From 



-dv 



[4m 2 (1 + u 2 ) + v 2 } 7 ' 2 v / T^ 2 15 (1 + 2m 2 ) 6 uVTT 

u 

we deduce that the original integral C 5 /640 is equal to 

-2 + 4 ((1 + 2u 2 f - l) E{iu) 



15tt 2 



/ 



-u du 



15tt 2 



15tt 2 



(1 + 2m 2 ) 5 

1 - 2 ((1 + 2u 2 f - l) +(1 + 2m 2 )] E(iu)K(iu) 



(l + 2u 2 )" 



u du 



00 

I 



-1 + ((1 + 2u 2 f - l) - (1 + 2m 2 )] K (iuf 



[l + 2u 2 )- 



u du 



which is (2/5)(C 4 /256). 



A more complicated but similar argument yields that C 6 = Ci- Applying the 
same procedure for n = 3, however, is fruitless. The hard integral in this special case 
becomes 

tt/2 tt/2 

J J sin (pyj sin 2 9 sin 2 p + cos 2 p — sin p dip d6 



tt/2 tt/2 

= J J cos p sin 2 (pyj sin 2 9 tan 2 ip + 1 — d9 dip 


tt/2 

= J cos p sin 2 ip E (i tan p) — dp 


00 

1 f t 2 E(it) 



4ttJ (l + t 2 ) 5/2dt 96 ^ 3 
v ' 

which bears no obvious resemblance to (8/3)(C 4 /256). Therefore C3 = C4 remains 
unproven, implying that ( 4 = 3n 3E2 (— |, |, |; 1, 2; l) is not yet a theorem. 
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